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ABSTRACT 

Cassiopeia A was observed using the Low-Band Antennas of the LOw Frequency ARray (LOFAR) with high spectral resolution. This allowed a 
search for radio recombination lines (RRLs) along the line-of-sight to this source. Five carbon a RRLs were detected in absorption between 40 and 

50 MHz with a signal-to-noise ratio of > 5 from two independent LOFAR datasets. The derived line velocities (v LSR 50 km s~') and integrated 

optical depths (~ 13 s ) of the RRLs in our spectra, extracted over the whole supernova remnant, are consistent within each LOFAR dataset 
and with those previously reported. For the first time, we are able to extract spectra against the brightest hotspot of the remnant at frequencies 
below 330 MHz. These spectra show significantly higher (15-80 percent) integrated optical depths, indicating that there is small-scale angular 
structure on the order of ~ 1 pc in the absorbing gas distribution over the face of the remnant. We also place an upper limit of 3 x 10~ 4 on the peak 
optical depths of hydrogen and helium RRLs. These results demonstrate that LOFAR has the desired spectral stability and sensitivity to study faint 
recombination lines in the decameter band. 

Key words. ISM: clouds - radio lines : ISM - ISM: individual objects: Cassiopeia A 



1. Introduction 



Radio recombination lines (RRLs) arise from ions recombining 
with electrons in diffuse, partially ionised gas. At high princi- 
pal quantum numbers (i.e. n ^ 200 at frequencies % 1 GHz), 
RRLs are an important probe of the temperature a nd density of 
the c o ol (T £ 100 K) int erstellar medium (e.g., iPavne et all 
[l989t iKonovalenkd fl984h . RRL measurements have distinct 
advantages over other tracers of the ionised gas; the extinc- 
tion caused by dust and the contamination from other sources 
are negligible. However, there have been only a limited num- 
ber of surveys that have studied RRL s from carbon at low 
radio frequencies (^ 330 MHz; e .g. Erickson et al. . [1995; 



Kantharia & Anantharamaiah, 2001; Roshi et al., 2002). Also 



the low angular resolution of these early surveys, from several 
degrees to around 6 arcmin, did not allow the location of the 
line-emitting gas to be determined, and also limited their capa- 
bility to detect faint RRL sources due to beam dilution effects. 



The LOw Frequency ARray (LOFAR; IStappers et aT1l2011l: 
Ivan Haarlem et al.ll2013l) is a new-generation radio telescope, 
which given its unique design, can be used to make sensi- 
tive surveys for RRLs with high spectral and spatial resolution. 
LOFAR comprises thousands of antennas arranged in (currently) 
41 stations spread across the Netherlands and other European 
countries. By design, LOFAR provides excellent baseline cov- 
erage, particularly in the inner 3 km core-region, where over 
half of the collecting area is contained. The low band antennas 
(LBAs) of LOFAR, although optimised for frequencies above 
30 MHz, cover the frequency range from 10 to 90 MHz. The 
high band antennas (HBAs) cover the frequency range from 110 
to 250 MHz. This large observable frequency window makes 
LOFAR an ideal instrument for carrying out studies of the ISM 
with low-frequency RRLs. 

In this Letter, we present the first spectral line measurements 
made using the LBA system of LOFAR. We target carbon RRLs 
in absorption towards the supernova remnant Cassiopeia A (Cas 
A), which is one of the brightest continuum sources in the sky 
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and has a sightline that is already kno wn to display strong car- 
bon RRLs in low frequency spectra (see lStepkin et al.ll2007l and 
references therein). In Sect. [2] we present our LOFAR commis- 
sioning observations and data analysis steps. Our carbon RRL 
detections are presented in Sect. [3] and we discuss these detec- 
tions and present our conclusions in Sect. |4] 

2. Observations & Data Reduction 

The data sets for our RRL study are taken from two LOFAR 
commissioning observations of Cas A in the 30 to 90 MHz 
(LB A) band. The first data set (Data-1) was taken on 201 1 May 
4 using the 16 core stations and 7 remote stations that were con- 
structed at the time, giving baselines up to about 24 km. These 
data were taken with 1-s visibility integrations and 244 separate 
frequency subbands, each with a bandwidth of 195 kHz and 256 
spectral channels. The second data set (Data-2) was observed on 
2011 October 15 using 24 core stations and 8 remote stations, 
and had a maximum baseline length of about 83 km. Data-2 was 
taken in a multi-beam mode that is used for standard continuum 
observations, with one beam on Cas A and a second beam on 
Cygnus A. This resulted in only half of the available bandwidth 
being used for the observation of Cas A. Also, a coarser spec- 
tral resolution of 64 spectral channels within each of the 122 
subbands of 195 kHz bandwidth was used. The visibility inte- 
gration time was again 1 s. After the first step of interfe rence 
removal using the AOFlagger routine (Offringa e t al.Ll2010l) . the 
data were averaged in time and/or frequency as required. 

In total, we produced 3 data sets for our work. First, to 
make a continuum image of Cas A, we produced an interme- 
diate data set from Data-1 that had a coarse spectral resolu- 
tion of ~ 50 kHz channeL 1 in order to speed up the calibration 
process later. Next, we created two RRL data sets. For Data- 
1, we included only those baselines between the separate core- 
stations (< 2.6 km) and retained the full spectral resolution of 
762 Hz channeL 1 . Data-2 was processed to retain a spectral res- 
olution of 3.1 kHz channeL 1 and to also include longer baselines 
(< 24 km) so that spectra could be extracted from individual re- 
gions across the supernova remnant. For these two RRL data 
sets, we decided to concentrate our search to those frequencies 
between 40 to 50 MHz. This was because the RRL absorption 
was expected to be stronger at the lower end of the LBA spec- 
trum due to the steep broadband spectral energy distribution of 
Cas A. The lower frequency end of the search window was cho- 
sen becau se the level of terre strial interference increases below 
40 MHz (Offringa et al., 2013). The resulting time and frequency 
resolution for each data set after these preprocessing steps are 
summarized in Table Q] 

We initially calibrated our continu um data using th e 
BlackBoard Self-calibration package (bbs; Pand ev et all [2009), 
with a relatively simple initial model, containing about a dozen 
point sources and one shapelet mo del. We then used the sage- 
cal package dKazemi et all 1201 lh to solve for the direction- 
dependent amplitude and phase variations in the data and to 
update the shapelet model for Cas A. To mitigate artefacts 
due to Cygnus A, present about 30 degrees away, we simul- 
taneously solved for direction-dependent gains towards Cas A 
and Cygnus A, and subtracted the latter in the visibility data 
using a shapelet model. We then imaged the source using 
the casa imager (http://casa.nrao.edu/) and constructed 
a composite sky model from the im age by fi t ting f or point 
and shapelet sources (see e.g., lYatawatta et all 1201 lh . We it- 
erated over this scheme (bbs calibration — » imaging — > sky- 
model — > calibration) to obtain high-resolution continuum im- 
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Fig. 1. LOFAR image of Cas A at 52 MHz (single subband of 
195 kHz bandwidth). The image was made using Briggs weight- 
ing (robust = 0) and was restored using a beam size of 40 arcsec. 
The hotspot region used to extract spectra from Data-2 is marked 
with the green square. 



ages and improved source models. In Fig.Q] we present a 40 arc- 
sec resolution continuum image of the supernova remnant at 52 
MHz (single subband). Given the excellent spatial resolution of 
LOFAR, a number of complex fea t ures a re clearly resolved (c.f., 
Kantharia et al I ll^inn7etdll200l . Further details of the 
calibration strategy and a full discussion of the continuum imag- 
ing of Cas A, and the wider surrounding field, will be presented 
in a companion paper (Yatawatta et al., in prep). 

To detect RRLs along the line-of-sight to Cas A, we cali- 
brated the RRL visibility data sets (see Table [TJ with bbs us- 
ing our final sky model that was obtained from our analysis of 
the continuum data. Multi-channel maps were constructed using 
the clean algorithm within casa, and image cubes from the in- 
dividual subbands were analysed separately. Spectra were then 
extracted over the whole remnant and over the brightest hotspot 
(see Fig.[T]i. The resultant spectra displayed a few artefacts, such 
as breaks at the subband boundaries and bad channels on the 
edge of each subband; typically around 5 per cent of the chan- 
nels per subband. These bad channels were removed during fur- 
ther analysis. Doppler correction terms are currently not incor- 
porated in the LOFAR imaging pipeline, hence the spectra were 
Doppler corrected after extraction. This leads to an uncertainty 
of S 1 kms 1 in the determination of the central velocities and 
line widths. This error is negligible compared to the velocity res- 
olution of our data. 



3. Results 

The subset of Data-1 that was used for the RRL analysis 
only contained the short baselines between the core-stations 
(^ 2.6 km), for which Cas A is practically unresolved. Therefore, 
spectra were obtained over the whole remnant with a spectral 
resolution of 762 Hz channeL 1 , which is equivalent to a veloc- 
ity resolution of ~ 5 km s _1 at 45 MHz. We report the detection 
(at > 16cr level) of four carbon a RRLs (C518a, C538a, C541a 
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Table 1. Summary of the LOFAR data sets used for our analysis: raw visibilities from Data-1 were processed twice to create two 
separate intermediate data sets for i) continuum imaging and ii) the RRL detections (see Section|2|. The velocity resolutions were 
computed at 45 MHz. 



Data set 


LOFAR id 


Start date 


Time 


Duration 




Processed resolution 








(UT) 


(h) 


Int. time (s) 


Channel width (kHz) Channel width (km s _1 ) 


Data-1 (continuum) 


L25937 


2011 May 4 


23:00:01.5 


18 


5 


48.768 325 


Data-1 (RRL) 


L25937 


2011 May 4 


23:00:01.5 


18 


5 


0.762 5.1 


Data-2 (RRL) 


L31848 


2011 Oct 15 


12:00:02.5 


15 


5 


3.1 21 



0.002 



0.000 



-0.002 



-0.004 



Cos A: 0(538) COMPARISON 



Cos A: WHOLE SNR vs. HOT SPOT 
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Fig. 2. A representative detection of a carbon a RRL towards Cas 
A with LOFAR. We show spectra of the RRL C538a that was 
measured using Data-1 (black solid) and Data-2 (blue dashed). 
For comparison, we have also smoothed the Data-1 spectrum 
to the same spectral resolution as Data-2 (red dashed). These 
spectra demonstrate the spectral stability of the LOFAR system 
for making spectral line measurements. 
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Fig. 3. The spectrum of the RRL C541a measured over the 
whole remnant (red) and only over the hotspot region (blue) us- 
ing Data-2. These data demonstrate that the integrated optical 
depth over the hotspot region is higher than across the remnant, 
showing that there is likely structure in the absorbing gas on 
scales of ~ 1 pc. 



and C543ff) towards Cas A in Data-1 between 40 to 50 MHz. 
In Fig. [2] we present a representative spectrum of a RRL detec- 
tion, with respect to the local standard of rest (LSR). For each 
detected RRL, we measured the line widths and integrated opti- 
cal depths (J t v dv) by fitting Gaussian components to the spec- 
tra, the results of which are presented in Table |2] Only single 
Gaussian-line profiles were used because our spectra do not cur- 
rently have a sufficient signal-to-noise ratio to reliably measure 
multiple components in the line structure. All of the detected 
RRLs have measured velocities at vlsr 48 + 1 kms -1 , con- 
sistent with the velocity of the gas in the Pe rseus arm, the likely 
absorber along the line-of-sight to Cas A QK onovalenko, 1984; 
lErshov et allll984ll987tlPavne et al.Lfl989Lll994 . These RRL 
detections from Data-1 represent the first spectral-line measure- 
ments made with LOFAR. 

The main goal of our study is to investigate whether there 
is any structure in the absorbing gas by detecting a variation 
in the integrated optical depths of the RRLs over the extent of 
the remnant. For this, LOFAR's unparalleled spatial resolution 
at frequencies below 100 MHz is required. We extracted spec- 
tra using Data-2 in two ways; one towards the compact hotspot 
component and the other over the entire area of the remnant. 
We detected five carbon a RRLs (at > 5<x level) in both our 
integrated spectra and the spectra extracted over the brightest 
hotspot component in the remnant. As these lines were only 
marginally resolved, we estimated their integrated optical depths 
by summing the measured optical depth over two channels. The 



results of these measurements are also given in Table|2] In Fig. [3] 
we show a representative detection of a carbon a RRL (C541a) 
again st the whole remn ant and the hotspot. Earlier studies (see 

e.g. JPavne et all[l9 89) reported two components at vlsr 50 

and ~ -40 km s _1 , where the former dominates by a factor 2-3. 
Due to insufficient signal-to-noise ratio we can not currently sep- 
arate these components, but the observed line-profile asymmetry 
in our spectra (see Fig. s [2] and [3} is consistent with two compo- 
nents. 

We also carried out a search for hydrogen and helium RRLs 
in the individual subbands of our data sets (without any fold- 
ing in frequency), but found no significant detection. The spec- 
tral noise per channel in our most sensitive data set (Data-1) 
is ~ 10~ 4 , so we derive an upper limit (3cr) on the peak op- 
tical depth of hydrogen and helium RRLs of 3 x 10~ 4 , or an 
equiva lent integrated optical depth of 0.9 s _1 . Following IShaved 
(1976), we derive an upper limit of 4 x 10~ I8 s~' for the in- 
terstellar ionisation rate of hydrogen. Our limit is significantly 
lower than the previously reported values despite t he uncertain- 
ties in determining various param eters (see e.g., iPavne et all 
1994: lGordon & SorochenkoTl2009L §3.3.6), indicating that cos- 
mic rays can not ionise the intercloud gas to the observed degree. 

4. Discussion & Conclusions 

In total, we have detected five carbon a RRLs towards Cas A 
from our two LOFAR commissioning data sets. In Fig. |4] we plot 
the integrated optical depths of these RRLs as measured from 
each LOFAR data set, and also the integrated optical depths for 
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Table 2. The derived line parameters for the RRLs measured 
from Data-1 and Data-2. For the case of Data-1, the spectral res- 
olution was sufficient to resolve the lines and the fitted Gaussian 
full width at half maximum (FWHM) is reported. For Data-2, 
the RRLs were not fully resolved and so the integrated optical 
depths were obtained over two spectral channel widths. The inte- 
grated optical depths have been calculated over the full remnant 
(SNR) for datasets Data-1 and Data-2, and around the region of 
the bright hotspot for Data-2 (see Fig. [TJ. 



30 



25 



Line 


Rest 


Center 


FWHM 


Integrated optical depth 




frequency 


velocity 


(kms- 1 ) 




(s- 1 ) 




(MHz) 


(km s _1 ) 


Data-1 


Data-1 


Data-2 






Data-1 




SNR 


SNR Hotspot 


C548a 


39.87 








12 ±2 14 ±2 


C543a 


40.98 


-48 ±1 


19 ±2 


12.8 ±0.5 


15 ±2 21+2 


C541a 


41.44 


-48 ± 1 


19±2 


13.3 ±0.5 


15 ±2 23 ±2 


C538ar 


42.13 


-48 ± 1 


20 ±3 


14.6 ±0.9 


13 ±2 24 ±2 


C518ff 


47.20 


-48 ± 1 


18 ± 2 


12.8 ±0.9 


11 ±2 20 ±2 



those RRLs that were previously measured towards Cas A with 
other telescopes (lErshov et all Il984l 119871: iPavne et all 1 19891 
11994 . as a function of the principal quantum number. We see 
that the integrated optical depths, as measured over the whole 
remnant, are consistent between Data-1 and Data-2 for the four 
RRLs that are common to both data sets. This demonstrates the 
spectral stability of our LOFAR commissioning data sets. We 
also find consistency for the RRL C^SSa, where we have both 
LOFAR and previously reported measurements in the literature. 
Although we currently do not have measurements for the other 
two previously detected RRLs (C502<? and C552a), the reported 
integrated optical depths of these lines are generally consistent 
with the RRLs that we have detected in the LOFAR data sets. 

Multiple carbon RRLs are usually modelled together in order 
to derive the physical conditions of the line-absorbing gas, since 
the electron levels are in non-local thermal equilibrium (non- 
LTE) conditions. Crucial inputs to these models include accu- 
rate line width measurements from high-spectral resolution data. 
Kantharia et al. (1998) show that previous low-frequency mea- 
surements of carbon a RRLs can be modelled as originating in 
cold (T e ~ 75 K), low density (n e ~ 0.02 cirT 3 ) clouds that are 
located approximately 1 15 pc away from Cas A. Higher density 
models (T e ~ 75 K, n e ~ 0.1 cm 4 ) are also able to fit the data, 
but these are less appealing due to their thermal pressures being 
very high. Although our current RRL detections are consistent 
with these models, given the agreement with the previous mea- 
surements, future observations of RRLs over the whole LOFAR 
observing window, and hence over a large and densely sampled 
range of principal quantum number, will differentiate between 
these competing models. The results from such a detailed anal- 
ysis using new observations will be reported in a forthcoming 
paper (Oonk et al., in prep). 

We also find from Fig. [4] that the integrated optical depths 
towards the hotspot are larger by up to almost a factor of 2, 
when compared to the integrated optical depths measured over 
the whole remnant. This shows that there is small-scale angu- 
lar structures and/or the geometry of the RRL-absorbing gas is 
such that there is var ying levels of absorptio n over the face of 
the remnant (see, e.g. Kanth aria et all Il998l) . This corresponds 
to structure in the RRL-absorbing gas on the order of ~ 1 pc. To 
model this structure, we plan to make new LOFAR observations 
with a higher spectral and spatial resolution in the future. 

Our results from the commissioning data sets demonstrate 
that LOFAR can be used to carry out spectral line observations 



9- 20 



• Whole SNR (Literature) 

• Whole SNR (Data-1) 

• Whole SNR (Data-2) 
» Hotspot (Data-2) 



500 510 520 530 540 
Quantum number n 



550 



Fig. 4. The integrated optical depths as a function of the prin- 
cipal quantum number for both data sets (Data-1 and Data-2) 
and for those previously reported in t he literature (lErshov et al. , 
1 1984 1 1 9871 IPavne et all 1 1 989L 1 1 994 . We find good agreement 
between the measurements made over the whole remnant from 
the two LOFAR data sets and with those in the literature. We also 
find that the integrated optical depths of the RRLs are higher 
over the bright hotspot component with respect to the mea- 
surements over the whole remnant, implying that there is likely 
small-scale structure in the absorbing gas distribution. 



at low radio frequencies. For the case of RRL studies, LOFAR 
will be an important tool in the future for two main reasons. First, 
given the up to arcsecond-scale resolution that can be obtained 
with LOFAR, beam dilution effects will not affect our detection 
sensitivity. Second, due to the large observable frequency range 
of LOFAR, from 10 to 250 MHz, hundreds of low-frequency 
RRLs will be detectable for the first time. Furthermore, due to 
the large bandwidth that can be observed during each observa- 
tion, many of these RRLs will be detected in a single observa- 
tion, making wide-field surveys for RRLs possible for the first 
time. The results reported here form part of the LOFAR Surveys 
key science project. In the near future, we plan to observe various 
bright sources (Galactic and extragalactic) with compact mor- 
phologies in the LOFAR LBA and HBA bands. Later studies 
will focus on RRLs towards regions of diffuse emission in the 
Galactic plane and fields with complex emission features. 
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